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ABSTRACT 

The normalisation of the matter power spectrum, Ojj* is an essential ingredient to predict the 
phenomenology of the low redshift universe. It has been measured using several methods, 
such as X-ray cluster counts, weak lensing and the cosmic microwave background, which 
have yielded values ranging from 0.7 to 1 .0. While these differences could be due to system- 
atic effects, they could also be due to physics beyond the standard ACDM model. An obvious 
possibility is the presence of non-Gaussian initial fluctuations in the density field. To study 
the impact of non-Gaussianity on each of these methods, we use a generalised halo model to 
compute cluster counts and the non-linear power spectrum for non-Gaussian models. Assum- 
ing scale invariance, the upper-limits on non-Gaussianity from the WMAP CMB experiment 
correspond to roughly a 4% shift in as as measured from cluster counts and about 2% shift 
through weak lensing. This is not enough to account for the current internal and mutual dis- 
crepancies between the different methods, unless non-Gaussianity is strongly scale dependent. 
A comparison between future X-ray surveys with a two fold improvement in cluster mass 
calibration and future cosmic shear surveys with 400 deg 2 will be required to constrain non- 
Gaussianity on small scales with a precision matching that of the current CMB constraints 
on larger scales. Our results argue for the presence of systematics in the current cluster and 
cosmic shear surveys, or to non-standard physics other than non-Gaussianity. 
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1 INTRODUCTION 

In the inflation-driven ACDM paradigm, the initial fluctuations 
are assumed to obey Gaussian statistics. Current observations 
of the cosmic microwave background (CMB) of galaxy clus- 
tering and of galaxy clu s ters are consistent w ith this hypothe- 
sis (eg. IVerde etaimOOlt iGaztanaga and Wagdl2003l) . However, 
models with non-Gaussian initial conditions ha ve been pro posed 
and a re not yet ruled out by observations (see lAvila-Reese et alJ 
120031 and reference therein). These models include ones based 

on in flation theor i es wi t h non-linearit ies of a single scalar 

field jMartin et alJ l200Ct iGangui et alJ l2002h or with multi- 
fields (Peebles 1999; Antoniadis et al. 1997; Linde and Mukhanov 
Il997t Bernardeau and Uzan 2002), or on cosmological defects 
iTurok and Spergelll990l;lAvelino et alll998l) . 

In this paper, we study how non-Gaussian initial conditions 
can affect the determination of C7g, the normalisation of the mat- 
ter power spectrum on 8 h Mpc scales. This normalisation is 
an essential ingredient to predict the phonemenology of the low 
redshift universe and has been measured using several methods. 
The recent CMB measurements with WMAP I Sperge l et all2003l) 
yield a constraint on the normalisation of o"g = 0.9 ±0.1 (68%CL). 
The abundance of galaxy clusters depends stron gly on oft and 
yields values in the range of <7g = 0.7-0.8 (eg. IPierpaoli et all 
2003) for £2,,, ~ 0.3. Weak lensing by large-scale structure, or 'cos- 



mic shear', provides a direct measurement of mass fluctuations 
and has recently been detected a nd measured (see R efregierl2003l 
Van Waerbeke and Mellier 2003, for recent reviews). The different 
cosmic shear surveys yield values of C7g in the range of 0.7 — 0.9 for 
the same value of £l m . The marginal discrepancies within and be- 
tween these different techniques may be due to residual systemat- 
ics, such as uncertainties in the mass-temperature of X-ray clusters, 
or the calibration of the shear in cosmic shear surveys. They could 
also be due to physics beyond the standard ACDM model, which 
has different effects on the various methods and surveys. In partic- 
ular, non-Gaussianity appears as a prime candidate to explain dis- 
crepancies between the different methods. Cluster abundance and 
cosmic shear indeed both probe non-linear structures at low red- 
shifts, but are sensitive to matter fluctuations on different scales. 
Non-Gaussianity may thus affect these two methods differently, 
while leaving the CMB power spectrum essentially unaffected. Re- 
ciprocally, it is interesting to establish whether the comparison be- 
tween the different methods with present and future surveys may 
constrain non-Gaussianity. 

To study the impact of non-Gaussianity on measure ments of 
(7s, we present a generalisation of the halo model 1 Ma and Frv 
2000 ; Seliak 2000). Building on the work of Robinson and Baker 
1 200C) for the mass function. lMa and F r^j2000 j) and lKwama et al. 
1 1999) for the halo bias] and iNavarro et alJ (l99r3) and 
Avila- Reese et alJl2003l) for the halo profiles, we compute the non- 
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linear power spectrum for arbitrary non-Gaussian models. This al- 
lows us to compute the cluster temperature function and the cosmic 
shear statistics for these models. We then study whether a discrep- 
ancy in the determination of tfg from the different observational 
methods can be explained by non-Gaussianity. We also discuss 
how future surveys can be compared to set constraints on non- 
Gaussianity. The paper is organised as follows. In Sj2| we present 
the different non-Gaussian models and discuss how they are con- 
strainted with current observations. In J3| we describe the non- 
Gaussian halo model. In i|4| we study the effect of non-Gaussianity 
on the determination of afc with the different methods. Our conclu- 
sions are summarised in $5\ 



2 NON-GAUSSIAN MODELS 

The statistics of the post-recombination linear density field p are 
usually described by first forming the density contrast S = (p — 
p)/p with respect to the average matter density p = (p) in the 
universe. We consider the probability distribution function (PDF) 
p(8,M) of the density contrast S smoothed with a spherical top- 
hat of radius R corresponding to an enclosed mass 
It is defined so that 



JdS p(8,M)=l, fd8p{8,M)8 = (8}=0. 
The variance of the distribution is given by 
/ dS p{8,M)8 2 = (8 2 ) = a 2 (M) 



(1) 



(2) 



and can be computed from the linear matter power spectrum f\i„ (k) 
using 



C7» = j 



d s k 



P lin (k)W 2 (k,M), 



(3) 



where W(k,M) = 3(s'mkR/kR - coskR)/(kR) 2 is the top-hat win- 
dow function. The lin ear power spectrum is calc ulated using the 
transfer function from lBond and Efstathiodll984h . 

To quantify the level of non-Gaussianity of the distribution, 
iRobinson et alj |200(J) have considered the excess probability a in 
the 3cr tail over that of a Gaussian PDF, i.e. 



a(M) = ■ 



(4) 



Another way of quantifying non-Gaussianity is the skewness 

fl 3 {M) = {8 3 } = f d8 p(8,M)8 3 (5) 

To quantify non-Gaussianity from CMB measurements, several au- 
thors (eg. Komatsu et al. 2003; Verde et al. 2001) have considered 
the non-linear coupling fy^ of the primordial (pre-recombination) 
potential field <5 to the Gaussian field <5g, defined as 

* = *g+/a<z.[*g-<* 2 >]- (6) 

This parametrisation is motivated by a class of inflation models. 
A conversion between fyi and the skewness 113 of the (post- 
combinatio n) density fiel d can be d erived using the results of 
iMatarrese et alj l2000l) andlVerde et alJfeOOlh. Current constraints 
from galaxy clusters lAvila-Reese et alJ 120031 : IRobinson et alJ 
2000) yield a < 4 at 95%C L. As we will see in CTl the recent 
WMA P CMB measurement (Komatsu et al. 2003; Robins orTet alJ 
2000) translates into a < 1.6 at 95%CL for M ~ 10 14 M R . 



a=1.00 
a = 5.02 
« = 9.61 
a=12.7 
a=13.8 




Figure 1. Log-normal distribution for several values of the non-Gaussian 
parameter a. The Gaussian case corresponds to a = 1, which is reached as 
J3 -> 0. For a = 5.02, a - 9.61, a - 12.7, and a = 13.8, the corresponding 
p values are 0.2, 0.4, 0.6 and 0.8, respectively. For PDFs with a variance of 
1 these correspond to ^3= 0.61, 1.32, 2.26 and 3.68. 



As a special case, we consider scale free models, where the 
shape of the distribution is independent of the scale R. In this case, 
the PDF can be written as 



p(8,M) = o{M)- l g ^cr(M) 5 
where g is a function which obeys 



L 



dxg(x) = l, / dxxg(x)=0, / dxxg(x) — l. 



[ 



Jo 
so as to verify equations J1I2L For the Gaussian case, 



g G (x) = (2n) *e 



(7) 



(8) 



(9) 



As an example of a scale invariant distribution, we consider 
the log-normal distribution which provides an adequate description 
of a range of non-Gaussian mode ls based on non-linear or m ulti- 
field inflation and cosmic strings iRobinson and Baker 200fJ). For 
this distribution, 



1 



8W = 



where 



2x(x + x )P 



exp 



In 



2 ( x+xp 
7 



2p 2 



?p 2 



(10) 



(ii) 



so as to enforce equation {S}. The parameter /3 controls the non- 
Gaussianity of the PDF and tends towards zero as the distribution 
becomes Gaussian. FigureQshows this distribution for several val- 
ues of j8 and corresponding values of a. 



3 GENERALISED HALO MODEL 

A number of tools exist for modelling the growth of structure since 
inflation. On the largest scales, fluctuations in density are suffi- 
ciently small that their growth can be accurately modelled using 
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linear theory. On small scales, however, the complex nonlinear in- 
teractions of gravitational in fall can, at present, only truly be mod- 
elled by using N-body simulation. These numerical simulations 
are however CPU intensive, making the exploration of parameter 
space prohibitively slow. To circumvent this, semi-analytic tools 
have been developed to provide fast and accurate fits to numerical 

simulations. 

( )ne such tool is the halo model jMa and FrvlbOOCt ISeliakl 
200G), which can be used to calculate the full non-linear power 
spectrum for a given cosmological model and has been shown to 
be in good agreement with high resolution numerical simulations. 
The halo model breaks the power spectrum into two parts: one that 
dominates on large scales and is determined by the clustering of 
halos of different masses, and a second term that is dominant on 
small scales and is determined by the density profile of the ha- 
los. The halo model thus requires three pieces of information: the 
mass-function of halos, the bias of halos compared to the back- 
ground density field and the halo d ensit y profile. The halo model 
as described bv iMa arid Fnl fcOOd) and lSeliald fcOOd) calculates 
these three ingredients by assuming that the density perturbations 
are Gaussian. We show how this can be extended to include models 
that have non-Gaussian fluctuations. 



3.1 Mass Function 

Robinson and Baker (2000) have proposed a generalisation of the 



Press and Schechter 1 1974) formalism which provides a good fit to 
N-body simulations for a class of non-Gaussian models. In their 
formalism, the fraction of objects of a given mass M or larger that 
collapse is calculated by integrating the PDF of the density fluctu- 
ations as 



F, 



F(>M) =A / P(S,M)d8, 



(12) 



where A is a normalisation factor and 8 C is a critical over density 
w hich can be com puted for a rbitrary cosmol ogies using the results 
o flLacev and ColeN 19931) and lEke et alii 199a) . The mass function, 
defined as the number density of halos per unit mass, is then given 
by 



dn p dv 
~dM ~~ MdM 
where 



/(v), 



/(v) 



dF 
dv 



and 



8 C 
a(M)' 



(13) 



(14) 



(15) 



The normalisation A is chosen to ensure that all the matter is ac- 
counted for, so that 



[p)i 



dn 

, dMM 

<p> Jo dM 



dv /(V) 



(16) 



For a scale invariant PDF (Eq. 0), it is easy to show that 
equations 1 12t and i 141 reduce to 

f(v)=Ag{v), (17) 

with A = [Jn'dx g(x)]~ . This simple form reveals that the mass 
function at mass M provides a measure of the normalised PDF g(x) 
at x = v = 8 c /<j(M). This assumes of course that the generalised 
Press-Schechter formalism is valid. 
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Figure 2. Mass function at z = for the log-normal model with different 
values of the non-Gaussian parameter a. The Gaussian case corresponds to 
a = 1, and the other values of a correspond to /3 = 0.01, 0.05, 0.2, and 0.4. 
For PDF's with a variance of 1 these in turn correspond to ^3=0.03, 0.15, 
0.61 and 1.32 



For the Gaussian case (Eq. (5)), we obtain the usual result 



/g(v) = \/ -e 
71 



-V-/2 



with the usual normalisation factor A 
IPress and Schechtedil974T) . 

For the log-normal distribution (Eq. 1101 ). we obtain 



(18) 

2 prescribed by 



/(v) 



with 



A = 



27r(v + x )j3 



exp 



ln 2 ( 



2 I V+.fQ N 



2/32 



(19) 



(20) 



where erfc is the complementary error function. Figure [2] shows 
the resulting differential mass function dn/dM for this model at 
z = 0, for several values of the non-Gaussian parameter a. Non 
Gaussianity with a > 1 tends to increase the number of halos with 
M > 10 14 M S and, to a smaller degree, to deplete halos with masses 
smaller than this limit. The high sensitivity of high mass halos to 
non-Gaussianity is expected since these halos probe the positive 
tail of the distribution. Throughout this paper we assume a standard 
ACDM model, with Cl m = 0.3, Q A = 0.7, h = 0.7, n = 1 and Y = 
0.21. For this figure we adopted Og = 0.77, but we will choose other 
values for this parameter depending on the data set considered. 



3.2 Halo Bias 

The bias parameter in the halo model describes how halos cluster 
compared to the underlying linear density field. General ising the 
formalism of iMo and White] <1996l) . lKovama et alJ Jl999h showed 
that the bias parameter for a general scale invariant PDF is given 
by 



d\nf{v) 



din v 



(21) 
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M (M ) 

Figure 3. Effect of a on the halo bias. The solid line shows the Gaussian 
case. We see that as we increase the level of fluctuations in the 3d tail, 
massive objects (high v) become less clustered, while less massive objects 
(low v) become more clustered. 



Note that this expression satisfies the condition J dv f(v)b(v) = 
1, which is required to ensure that the lin ear power spe ctrum is 
recovered on large scales (see discussion in Seliak 2000). For the 
scale invariant case, /(v) is proportional to the PDF (see i]3.1> . The 
halo bias can therefore be thought of as a measure of the derivative 
of the PDF in this case. The measurement of halo bias, using the 
halo correlation function for instance, can thus be used to probe the 
PDF independently of the mass function. 

For the Gaussian case (Eq. 1181 ). this reduces to the usual Mo 
& White result 



Mv) = i- 



(22) 



More accurate bias functions (eg. I Jin J 19981) have been devised to 
provide better fits to N-body simulations, but they cannot be readily 
generalised to non-Gaussian models. Since we are only interested 
in the effect of non-Gaussianity relative to the Gaussian case, the 
above expression suffices for our purposes. 

For the log-normal distribution, the bias function is 



b(v) = l + 



XQ 



P 



■In 



V+x 



V 



(23) 



which does indeed go to Eq. 1221 as j8 tends towards zero. Figure 
|5|shows the bias as a function of mass for several values of a. For 
skew-positive PDFs (a > 1), objects with low (high) masses are 
more (less) clustered compared to the Gaussian case. 



3.3 Halo profile 

We model the density profile of the halos as 

p(r)=p s u(r/r s ), u{x) = x^{\+xf- i (24) 

where r s is a characteristic radius. The inner slope p must be 
less than 3 for convergenc e, a nd is equal to 1 in the case of a 
Navar ro, Frenk & White Jl99fl NFW) profile and 1.4 in the case 
of a iMoore et alJ J 19981) profile. We fix the normalisation p s so 



that the density contrast within the viral radius r v is S v = 200, 
i.e. so that M = 4n~p8 v rl/3 = 4n /q" dr r 2 p(r). This yields p s = 
p~8 r c i U~ l (c)/3 where we have defined the compactness parame- 
ter c = r v /r s , and 

U(c)= dxx 2 u(x) = 2 Fi(3~p,3-p-A-p;-c) (25) 

Jo 3 - p 

where 2 F\ is the hypergeometric function which can be evaluate 
numerically I Pr ess et alll996l) . In the case of the NFW profile (p = 
1), this function reduces to U(c) = ln(l +c) — c/(l +c). 

For the ACDM Gaussian model, a good fit to the non-linear 
power spectrum is given for an NF W profile (/X = 1) wi th the 
concentration parameter is chosen as ICooravlExxjllCoorav et ail 
2000) 



c(M,z) =c*(z) 



M 



M*{z) 



-Hz) 



where c*(M,z) ~ 10.3(1 



-0.3 



and h(z)~ 0.24(1 



-0.3 



(26) 



Here 



(z) is the non-linear mass s cale defined by v(M„,z) = 1. 
lAvila-Reese et alj (2003) have used numerical simulations to 
measure the halo profiles in non-Gaussian models. For a model 
with a ~ 3, they find an increase in the compactness parameter c 
of about 20% compared to the Gaussian case for halos with masses 
in the range 10 11 to 10 15 h M©, They also find the inner slope to 
be steeper (p ~ 1.6) compared to the Gaussian NFW case (p = 1). 

More simulations would be needed to derive new scaling rela- 
tions for halos in general non-Gaussian models. For our purposes, 
we use the use the NFW profile as a reference, but we also study the 
ef fect o f changin g the comp actness and inner slope to those found 



ettect ot changing tne compai 
bv lAvila-Reese et alJ (2003 ) . 

3.4 Non-linear Power Spectrum 
iMa and Frvl <200fJ) and lSeliakl i20o"i 



showed that the halo model 
provides an accurate prescription to compute the non-linear power 
spectrum for Gaussian models. In this formalism, the non-linear 
power spectrum P (k) at a given redshift is written as 



P(k)=P l (k)+P 2 (k), 
where the 1-halo term is 



f°° dn 

PAk) = dM 

v ; Jo dM 



p(k,M) 



and the 2-halo term is 
dn 

dM 

lo 



m) = 



dM —b(M) 
dM v ' 



p(k,M) 



PUk)- 



(27) 



(28) 



(29) 



In these expressions, p(k,M) is the radial Fourier transform of 
the density profile (Eq. 1241 ) of a halo of mass M. The de- 
tails of our impl ementation of the halo model can be found in 
iRefregier and Tevssierl l2002t) . 

Figure [4] shows the effect of non-Gaussian initial conditions 
on the non-linear power spectrum. Non-Gaussianity is seen to af- 
fect transition scales 0.1 < k < 3 h Mpc -1 between the linear and 
the non-linear regime. The linear (k < O.lh Mpc -1 ) and strongly 
non-linear (k > 3h Mpc~') are not significantly affected. Figure|5| 
shows the effect of varying the inner slope p and the concentra- 
tion parameter c as ind icated by the non-Gaussian simulations of 
lAvila-Reese et alJj2003h . Their suggested alteration of the concen- 
tration parameter by 20% (which we applied to all masses) tends to 
only slightly increase the power spectrum on small scales (k > 3h 
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Figure 4. Non-linear power spectrum at z = for several values of a. The 
a = 1 Gaussian model (solid line) is a standard ACDM model with Q. m = 
0.3, r = 0.21 and (Jg = 1.0. Its decomposition into the 1-halo and the 2-halo 
terms are shown as long dashed lines. 

Mpc -1 ). On the other hand, the alteration of the inner profile /i 
by 60% increases the power spectrum significantly on these scales. 
Note however, that the mass dependence of c(M) which we used as 
a reference (Eq. 1261 ) was chosen to reproduce the non-linear power 
spectrum in numerical simulations for an NFW profile (fl = 1). 
A selfconsistent fit to the non-linear power spectrum from non- 
Gaussian simulations would be needed to establish whether this 
effect is indeed so pronounced. In the following, we keep c and fl 
unchanged from the Gaussian NFW case, but remain cautious in 
the interpretation of our model for the power spectrum on small 
scales. 



4 POWER SPECTRUM NORMALISATION 
4.1 Cosmic Microwave Background 

Using the recent measurements from the WMAP mission, 
ISpergel et alj <2003l) found ct 8 = 0.9 ±0.1 (68%CL). When com- 
bined with other CMB missions, galaxy surveys and Lyman al- 
pha surveys, the authors find <7 8 = 0.84 ± 0.04 (68%CL). The 
CMB power spectrum provides a 2-point statistics of the matter 
fluctuations in the linear regime and is thus not affected by non- 
Gaussianity 1 This is also true when it is combined with measure- 
ment of the galaxy power spectrum on linear scale, while assuming 
that the galaxy bias is scale independent. As a result, the determi- 
nation of Og based on fits to the CMB power spectrum (with and 
without galaxy surveys) is insensitive to non-Gaussianity. 

The measurement of higher order statistics of CMB maps, 
however, does currently present tight constraints on the level 
of primordial non-G aussianity. Using the WMAP measurements, 
iKomatsu et aljl2003l) were able to place limits on the non-Gaussian 

1 Strictly speaking, non-Gaussianity may slightly bias the algorithms used 
to measure the CMB power spectrum, which often assume exact Gaussian 
statistics. In this paper, we ignore this small effect. 



Figure 5. Effect of the halo profile on the non-linear power spectrum. For 
a = 3, Avila-Reese et al. (2003) find an increase in the compactness param- 
eter, c„ g that is roughly 20% greater than that for the NFW profile, Cnfw- 
They also find the inner slope /i„ g ~ 1.6 to be more cuspy for this level of 
non-Gaussianity, than that for the standard NFW profile Unfw = 1. 

coupling defi ned in equation j(p of - 58 < /nl < 134 (95%CL). 
According to Matarr ese et alJ OOOOl figure 2), the upper limit 
corresponds to an upper limit on the skewness of the post- 
recombination density field (Eq. (5J) on scale of M = 10 14 M s of 
i$ 8.0 x 10~ 2 , for a ACDM model with ofc = !• For the log- 
normal distribution (Eq. 1101 ). this corresponds to a < 1.6 at the 
95% confidence level. Note that this conversion assumes that the 
scale dependence of the non-Gaussianity is that of models given by 
equation 

4.2 X-ray Cluster Counts 

Counting X-ray clusters is a direct probe of the mass function and 
provides a sensitive measure of o"g. The main difficulty is to es- 
timate the masses of X-ray selected clusters. The most reliable 
method to do so consists of measuring the X-ray temperature of 
the cluster and to use a mass-temperature (M-T) scaling relation. 
Using an M-T relation derived form hydrodynamical simulations, 
several groups have derived a value of <jg around 0.9 from mea- 
surem ents of the X-ray temperature function (see IPierpaoli et all 
2001 , and reference therein). More recently, the use of an observa- 
tio nal M-T r elation led to a revision of erg in the range 0.7-0.8 (see 
eg. Pierpao li et all2003llSeliakl200ll) . The uncertainty in the M-T 
relation is, at present, the largest source of uncertainty in cosmo- 
logical parameter estimates using this method. 

In order to study how the determination of c>8 using 
this method is affected by non-Gaussianity, we consider the 
compilation of clust er temperature measurements presented by 
Pierp aoli etaflfcOOlh . The catalogue consists of temperature mea- 
surements for 38 clusters taken using ROSAT and ASCA. The clus- 
ter temperatures range from kT = 3.5 to 12 keV, with those with 
temperatures greater than 7.5 keV forming a volume limited sam- 
ple. The resulting cumulative X-ray temperature function n(> kT) 
from this compil ation is show n in figure|6| To compute the cluster 
masses, IPierpaoli etidlfcOOll) have used the M-T relation, 
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Figure 6. Observed temperature function from the compilation of Pierpaoli 
et al. (2001). Also shown is the predicted temperature function for several 
values of the non-Gaussian parameter a and for o~8 = 0.77. 



Figure 7. Determination of o~g from a fit to the cluster temperature function 
as a function of the non-Gaussian parameter a. Several temperatures kT for 
the X-ray cluster counts are displayed. 



M(T,z) 



A C E' 



2S1 A 

^7 



(30) 



where T is the temperature in keV, and E = £X m (l + z) + &a + 
+ z) 2 . The parameter A c is the mean over density inside a 
vir ial radius and c a n be c omputed using the fitting fun ctions given 
bv IPierpaoli et alj 1200 ll) . In a later work. Ipieroaoli et"all 1 2003 ) 
fit the observed temperature function using this relation and find 
(7 8 = 0.77+°$ (68%CL). 

Figure|5|also shows the predicted temperature function for dif- 
ferent values of the non-Gaussian parameter a. For the masses cor- 
responding to kT > 1 keV (corresponding to M > 3 x 10 13 /i _1 M e ), 
a skew-positive PDF (a > 1) enhances the cluster counts (see Fig- 
ure the effect being more pronounced for larger temperatures. 
As a result, the determination of CTs from the X-ray cluster counts 
will be affected by non-Gaussianity. Figure shows the value of 
(78 derived fr om a fit to the temperature function as a function of 
a. Following Seliak 1 2001), we considered a fit to n(> kT) at a 
single temperature kT, which we choose to be 4, 7, and 10 keV, 
in turn. For kT ~ 7 keV and for the Gaussian case ( a = 1), we 
find ff a ~ 0.77 in agreement with lPierpaoli et all 12003) and lSeliakl 
1200 ll) . The value of a% required to fit the temperature function de- 
creases as a increases. The effect is more pronounced for larger 
cluster temperatures, as expected. 



4.3 Cosmic Shear 

Weak lensing by large scale structure, or cosmic shear, pro- 
vides a dir ect measure of the mass distrib u tion in the uni- 
verse (s ee [Bartetoiaimand Schneidei] 1200 it iRefresieJ 120031 : 
IVan W aerbeke and Mellier 2003j for recent reviews). This ef- 
fect is based on the shear that intervening large-scale struc- 
tures induce on the images of background galaxies. Several 
groups have now perfo r med accurate measurements of cosmic 
shear Ijarvis et alJl2003t lHamana et alj|2002t ferown et alJl2003t 
iHoekstra et all l2002t iRefresier et all l2002t Ivan Waerbeke et all 



l2002tlBacon et alJl2002l) . which have been used to set constraints 
on cosmological parameters. 

The 2-point cosmic shear st atistics can be desc ribed by the 
weak lensing power spectrum (see B acon et all 120001) for conven- 
tions) 

2 
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■X 



(31) 



where I is the multipole moment, % is the comoving radius, r is the 
comoving angular-diameter distance, a is the expansion parameter, 
and Ho and £l m are the present value of the Hubble constant and 
matter density parameter, respectively. It can be computed for non- 
Gaussian models by using the halo-model approximation for the 
non-linear power spectrum P(k,x) at a given comoving radius %, 
given in equation 1271 . The radial weight function is given by 



[ dx'n{x 

■17 



„ r( x)rjx'-X) 
r(x') ' 



(32) 



where n(x) is the red shift distribution of galaxies and is normalised 
asfdxn(x) = 1. 

Figure [8] shows the weak lensing power spectrum for sev- 
eral values of a. The galaxy redshift distributions was taken to be 
n(z) <x z 2 exp(— (z/zo) 1 ' 5 ) with zo = z m /1.41 and the median red- 
shift was chosen to z m = 1 as suitable for several cosmic shear sur- 
veys. The cosmological parameters are those for our fiducial model 
given in i!3.1l with o"8 = 1.0. The linear power spectrum is also plot- 
ted for comparison and reveals that non-linear corrections are very 
important on scales £ > 10 3 . Non-Gaussianity has the most pro- 
nounced effect on intermediate scales 10 2 < t < 10 4 , which is con- 
sistent with its effect on the 3D power spectrum (see figure|4j- For 
a > 1, non-Gaussianity tends to increase the lensing power spec- 
trum on these scales. 

In practice, it is often more convenient to consider the shear 
variance Oy = (|y| 2 ) in randomly place circular cells of radius 6. 
This is related to the shear power spectrum by 



diC e \W e [ 



(33) 
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Figure 8. Weak lensing power spectrum Q for several values of the non- 
Gaussian parameter a. The galaxy median redshift was assumed to be Zm — 
1. The Gaussian model (a = 1) is our fiducial ACDM model with o~g = 1. 



Figure 9. Shear variance <Jy(6 ) as a function of the circular cell ra- 
dius. Both recent measurements jBrow n et all2003: Refregi er et al]|2002l 
IVan Waerbeke et alj|2002t iBacon et al]|2002i) and the predictions for the 
different non-Gaussian models (with z m = 1 and °~8 = 1) are shown. 



where Wi = 2J\ (£6) /(£6) is the Fourier transform of the cell aper- 
ture. The shear variance of the different non-Gaussian models is 
shown in figure |9| Also shown are recent measurements of the 
shear variance from different surveys with median galaxy redshifts 
in the range 0-8 < z m < 1 iBrown et all2003l:|Refregier et all2002l 
IVan Waerbeke et alj|2002t IBacon et alj|2002h . The measurements 
are in rough agreement with each other and with the Gaussian 
ACDM case. Non-Gaussian models with a > 1 produce larger 
shear variances on scales 6 < 10', i.e. on scales to which most 
surveys are sensitive. This is not surprising since this limit corre- 
sponds to the transition between the linear (9 > 10') and non-linear 
(6 < 10') regimes for the shear statistics. 

We therefore expect non-Gaussianity to affect the determina- 
tion of rTg from current cosmic shear surveys. Figure fTol shows the 
value of <7g derived from the shear variance measured at different 
scales and with z m = 1, as a function of the non-Gaussian parameter 
a. We see that if the universe is non-Gaussian with, say, a ~ 3, the 
cosmic shear surveys (which assumed a=i) would have overesti- 
mated o% by about 5%. This shift is approximately unchanged if the 
galaxy redshift is set to z m =0.6, as appropriate for shallow er cos- 
mic shear surveys irarvis et all2003t hoekstra et alfcooj), which 
tend to find lower values of O r than d eeper ones iBrowri^Hil .2003; 
iRefreeier et all2002tlvan Waerbeke et all2002tlBacon et all2002t 
see however Hamana et al. 2002). The discrepancy between shal- 
lower and some of the deeper surveys can therefore not be ex- 
plained by non-Gaussianity. 

As discussed in 33.31 the effect of non-Gaussianity depends on 
our assumption regarding the halo profile. We find that the bias in 
fjg is twice as pro nounced if we use the modified profile suggested 
by lAvila-Reese et alJ <2003l) instead of the NFW profile. Further 
high-resolution numerical simulations would be needed to resolve 
this theoretical uncertainty. 

As expected, the effect of the non-Gaussian parameter, a, de- 
pends on the scale being studied. On large angular scales, 6 > 10', 
a has a very weak effect on the shear variance, which is as expected 
since these scales are still in the linear regime. The intermediate 
scales, 1' < 6 < 10', are the most affected by our non-Gaussian 
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Figure 10. Determination of o~8 from a fit to the cosmic shear variance 
Gy (9) as a function of the non-Gaussian parameter a. Several scales 9 for 
the cosmic shear variance are displayed. 



PDF, as can also be seen in figure fTol This is also consistent with 
the behaviour of the weak lensing power spectrum (figure[8j since 
these intermediate scales correspond to 10 2 < I < 10 3 . On scales 
smaller than 8 = 1', the variance from models with different a val- 
ues no longer diverge from each other. This makes measurements 
on small angular scales less sensitive to variations in the PDF. On 
even smaller scales, corresponding to I > 10 4 , the shift in the power 
spectrum changes sign, suggesting that on these scales the spread in 
shear variance of our models will decrease. Since lensing surveys 
do not probe these angular scales, we do not investigate this further. 
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Method 


Scale 


s 


°f 

o 


cluster 


4.0 keV 


0.032 


0.74 


counts 


7.0 keV 


0.043 


0.77 




10.0 keV 


0.052 


0.82 


weak 


0.5' 


0.014 


0.99 


lensing 


1.0' 


0.017 


0.99 




5.0' 


0.019 


0.99 




10.0' 


0.015 


0.99 




20.0' 


0.011 


0.99 



Table 1. Fit parameters for the error Ao"8 = (Jg — o~g ~ s(a — 1) made in 
the measurement of o~8 if non-Gaussianity is ignored. 

4.4 Combined Constraints 

As discussed in the previous sections, the determination of <7g from 
the CMB with and without galaxy surveys is insensitive to non- 
Gaussianity. On the other hand, cluster counts and cosmic shear 
are affected by non-Gaussianity, as can be seen in Figures [7] and 
1101 For a < 5, we find that the error in (7g made by ignoring non- 
Gaussianity can be approximated by 

AcT 8 = <T 8 G -cr 8 ~. s (a-l), (34) 

where o% is the value derived by (wrongly) assuming a Gaussian 
model (i.e. corresponding to a = 1 in these figures). The slope 
i' and intercept values 0% are listed in table Q for cluster counts 
and cosmic shear at different temperatures and angular scales. 
As is apparent from the figures and in the table, cluster counts 
are more sensitive (s larger by a factor of about 2) than cosmic 
shear to non-Gaussiani ty. However, as noted above, the use of the 
lAvila-Reese et alJ 120031) halo profile instead of the NFW profile 
would tend to increase the cosmic shear slope s by about a fac- 
tor of 2, which would then nearly match that for cluster abundance. 
While this needs to be established by further numerical simulations, 
we will use the NFW slopes s for cosmic shear as a reference but 
remember that they are lower limits. 

If we assume that the PDF shape is scale invariant over all 
scales, we have seen (see i|4. It that the recent WMAP skewness 
measurement imposes a constraint of a < 1.3 (68%CL). This limit 
corresponds to errors in <y$ of A<7s ~ 0.02 and 0.01, for clus- 
ter counts and cosmic shear respectively. These are significantly 
smaller to both the current uncertainties (Atfg ~ 0.05 and 0.1, 
respectively) and the discrepancy between the different methods 
and measurements (Acts ~ 0.3). Thus, with this assumption, non- 
Gaussianity can not explain the internal and external discrepancies 
between these methods. 

Note however that the above conclusion rests on an extrapo- 
lation from scales probed by WMAP (1 < 1000 corresponding to 
> 25 h~ l Mpc comoving at decoupling) down to scales relevant 
for clusters and cosmic shear, 0.1 to 5.0 comoving /i^Mpc. If wc 
drop the assumption that the PDF is invariant over this large range 
of scales, the WMAP measurement no longer constrains the PDF 
relevant for these two latter techniques. In this case, we can attempt 
to reconcile the cluster and cosmic shear normalisations using non- 
Gaussianity. We find that this would require jS ~ —0.4. This level 
of negatively skewed PDF does not have a corresponding a param- 
eter since the distribution does not have a 3fJ tail in this case. Such 
a negative skewness would be inconsistent with the observed shape 
of the temperature function (see figure[6) and of the shear variance 
(see figure^- 



Finally, we look at the level of accuracy future surveys need 
to reach in order to improve upon the current constraints on non- 
Gaussianity. Assuming scale invariance, we have shown that the 
current la constraints set by WMAP correspond to uncertainties 
of AfJg ~ 0.02 for cluster counts and AO's ~ 0.01 for cosmic shear. 
At present, cluster count measurements have Ict uncertainties of 
AcTg ~ 0.05 which are dominated by an uncertainty in the normali- 
sation r* of the M-T relation of about 10%. This means that if the 
error on fTg from clusters is reduced by a factor of about 2 then they 
will be comparable to that arising from our lack of knowledge of 
the statistics of the density field. This will be achieved if the un- 
certai nty in T* is red uced to about 4% in future surveys (see eg. 
IPierpaoli et ail koOS). Similarly, current lensing surveys of about 
10 de g 2 have 1<7 errors of Acts ~ 0.06 (eg. IVan Waerbeke et alJ 
2002 ) and scale as \/A where A is the survey area. Lensing surveys 
will thus have to increase in sensitivity by a factor of about 6, and 
thus have areas of the order of 400 deg 2 , before their measurement 
of tfg becomes sensitive to non-Gaussianity. Note, however, that it 
may be possible to probe primordial non-Gaussianity by studying 
higher order statistics of cosmic shear. 



5 CONCLUSIONS 

We have investigated the impact of non-Gaussian conditions on the 
determination of CTs via various observational methods. For this 
purpose, we have generalised the halo model to compute cluster 
statistics and the non-linear matter power spectrum. In passing, we 
noted that, for scale free non-Gaussian models, the halo mass func- 
tion and bias provide a measure of the linear density PDF and of its 
derivative, respectively. The measurement of both of these statis- 
tics, therefore, provides an independent measure of the PDF and a 
test of its scale invariance. 

The determination of CTs from the CMB power spectrum and 
from galaxy surveys in the linear regime is not sensitive to non- 
Gaussianity. On the other hand, the cluster X-ray temperature func- 
tion tends to be enhanced at large temperatures (T > 3 keV) for 
skew-positive (a > 1) non-Gaussianity. If the primordial fluctua- 
tions are non-Gaussian and this is ignored, the determination of 
CTs from current X-ray cluster surveys would be overestimated by 
Acts ~ 0.04(o: — 1). This depends weakly on the typical tempera- 
ture of the clusters. Cosmic shear statistics from current surveys 
probe the non-linear part of the matter power spectrum and are 
thus also sensitive to non-Gaussianity. Using the generalised halo 
model, we find that the matter power spectrum is enhanced by a 
factor of 2 on scales OA < k < 3 h Mpc -1 for skew-positive non- 
Gaussianity of a = 5. As a result, the cosmic shear power spectrum 
is also enhanced by a factor of 2 on scales 10 2 < i < 10 4 . The error 
in o"8 from current cosmic shear surveys if non-Gaussianity is ig- 
nored is Acts ~ 0.02(a — 1). This behaviour depends weakly on the 
angular scale 9 and on the galaxy redshift z m - We also find that the 
nonlinear power spectrum, and therefore cosmic shear surveys may 
be sensitive to the halo profile. This, in turn, may be sensitive to the 
primordial density field and needs to be st udied in m ore d etail v ia 
high-resolution N-body simulations (see lAvila-Reese et alj2 003). 

The strongest upper limits on primordial non-Gaussianity are 
provided by the recent WMAP measurements and correspond to 
a — 1 < 0.6 (95% CL). Within this limit and assuming scale invari- 
ance, non-Gaussianity cannot explain the differences between the 
CTs values derived from different cosmic shear surveys and X-ray 
cluster catalogs. Dropping the assumption of scale invariance from 
the scales probed by WMAP to those probed by cluster counts and 
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shear, non-Gaussianity would tend to decrease the derived values 
of (7g for both cosmic shear and cluster counts and thus would need 
to be very pronounced to explain the discrepancy. Specifically, a 
highly negatively skewed PDF with /3 —0.4 (which does not have 
a corresponding a) would be needed to resolve the discrepancy, 
but would not be compatible with the observed shape of the clus- 
ter temperature function and of the cosmic shear 2-point function. 
A comparison between future X-ray surveys with a two fold im- 
provement in cluster mass calibration and future cosmic shear sur- 
veys with 400 deg 2 will be required to constrain non-Gaussianity 
on small scales with a precision matching that of the current CMB 
constraints on larger scales. Our results suggest that the discrepan- 
cies are either due to systematics in one or several of the methods 
or to non-standard physics other than non-Gaussianity. 
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